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Structural  Mechanics;  Task  240104,  Vibration  Prediction  and  Control,  Measurement 
and  Analysis;  Work  Unit  24010432,  Large  Space  Structures  Technology  Program. 
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Section  I 


Introduction 

The  Structural  Vibration  and  Acoustics  Branch  has  initiated  a  program  to 
study  the  dynamics  and  control  of  Large  Space  Structures  (LSS).  This  Large 
Space  Structures  Technology  Program  (LSSTP)  is  intended  to  enable  the 
Flight  Dynamics  Laboratory  to  instrument,  test,  and  analyze  large  space 
structures  on  the  ground  in  order  to  predict  their  behavior  in  space.  Since 
the  testing  is  to  be  done  in  a  ground  based  laboratory,  i.e.  under  1-g 
acceleration,  the  experiments  must  be  designed  so  as  to  counteract  this 
gravitational  effect.  One  proposal  is  to  use  soft  suspension  systems.  Long 
cables  can  provide  pendulum  support  with  low  frequency  for  horizontal 
motion,  but,  for  a  pendulum,  the  restraint  in  the  vertical  direction  is  rigid. 

One  way  to  provide  soft  restraint  vertically,  while  maintaining  the  pendulum 
approach,  is  to  counter  balance  the  test  model. 


The  simplest  counter-balanced  suspension  system  is  the  classic  Atwood's 
machine. [1]  This  consists  of  two  masses  m,  and  m2  connected  by  an 

inextensible,  flexible  wire  of  negligible  mass,  draped  over  a  frictionless, 
massless  pulley,  which  in  turn  is  rigidly  suspended  from  an  overhead  support. 
The  motion  of  the  two  masses  is  assumed  to  be  constrained  to  the  vertical 
direction:  either  there  is  no  motion,  the  situation  that  occurs  when  the  two 
masses  are  equal  and  there  is  no  initial  velocity,  or,  as  one  mass  rises  the 
other  falls. 


In  previous  reports  the  Atwood's  configuration  was  studied  when 
the  one  mass,  mr  remained  constrained  to  move  vertically  but  the  other 

mass,  m2,  was  swung  as  a  pendulum.  The  results  described  in  these  two 

reports  may  be  summarized  as  follows:  after  deriving  the  equations  of 
motion,  these  were  programmed  and  run  on  a  digital  computer  using  a 
Runga-Kutta-Fehlberg  ^  routine.  The  first  numerical  experiments  assumed 
the  two  masses  m,  and  m2  were  equal.  As  one  normally  does  with  a 
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pendulum,  the  mass  m2  was  displaced  through  an  initial  angle  60  ,  released 
with  zero  initial  velocity,  0o(O)  =  0  and  the  subsequent  motion  studied. 
Several  runs  were  made  ,  taking  successive  values  for  60  of  0.01 , 0.10,  0.20, 
and  0.50  radians.  In  all  four  cases,  the  mass  m1  initially  dropped  ever  so 

slightly  and  then  rose  monotonically.  Obviously,  this  is  equivalent  to  a 
lengthening  of  the  pendulum  arm.  As  the  pendulum  arm  increased,  the 
amplitude  of  swing  decreased. 

While  not  concerned  per  se  with  counter-balanced  suspension 
systems,  material  of  collateral  interest  may  be  found  in  [5], [6]. 

These  two  references  treat  problems  whereby  a  swinging  mass  is 
being  hoisted  by  a  winch. 

It  is  intuitively  obvious  that  if  the  mass  m1  is  made  slightly  greater  than 
m2  that  this  monotonic  lengthening  of  the  pendulum  arm  with  time  would  be 
slowed  and  that  if  m1  were  made  much  greater  than  m2,  the  length  of  the 

pendulum  arm  would  actually  decrease.  Thus,  there  must  be  some  critical 
relation  between  and  m2  for  which  the  length  of  the  pendulum  arm  is 

either  constant  or  a  periodic  function  of  time.  Further  numerical 
experiments  [2]  indicated  that  this  relationship  is: 

m1=m2[l  .000+7.451x1  O'8  0o+O.25O  0Q2] 

It  was  found  using  a  least  square  fit  to  data  points  corresponding  to  0o=O.1 , 
0.2,  0.3,  0.4,  and  0.5  radians. 

In  this  report,  a  theoretical  derivation  of  the  relationship  between  m1  and 
m2  is  obtained  under  the  added  restriction,  however,  that  the  initial  angular 


displacement  0Q  is  small.  The  technique  used  is  classical  perturbation 
theory 


Section  II 


The  Mathematical  Statement 

The  geometry  of  the  configuration  studied  is  given  in  Figure  1 .  From  this 
geometry,  assuming  that  the  masses  of  the  pulley  and  the  wire  are  negligible 
and  that  the  radius  of  the  pulley  may  be  neglected  also,  the  differential 
equations  describing  the  motion  of  the  mass  m2  are: 

p  0  +  2pQ+  g  sin0  =0 

(m1  +  m2)  p  -m2p02  +  m1g  -  m2gcos0  =  0 

where  dots  indicate  differentiation  with  respect  to  the  time  t.  Equations  1 
and  2  may  be  simplified  slightly  if  we  observe  that  by  letting  p  =  pg  ,  both 
equations  contain  g  as  a  factor  and  we  may  divide  through  by  it.  This  is 
equivalent  to  choosing  units  so  that  g  *  1 .  Furthermore,  by  using  the  small 
angle  approximation,  these  equations  then  become: 

p  0  +  2p0  +  0  =  0 

(m1  +  m2)  p  -  m2p02  +(m1-m2)  +  !^m202  =  0 
Case  i)  Constant  p. 

We  prove  now  that  there  is  no  solution  of  Equations  3  and  4  for  which 
p  =pQ=constant.  Actually,  we  shall  show  that  assuming  that  p  =  constant 
leads  to  a  contradiction.  Using  this  assumption,  Equations  3  and  4  become 

p0  0  +  0=0 

m2pQ02-  (m1-m2)  -  v2m2Q?  =  0 
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Differentiate  Equation  4.1  with  respect  to  the  time  t  and  divide  through  by 
m2  to  get: 

2p0  00  -  00  =  0 

If  0  =  0,  then  Equation  3.1  implies  that  0=0  and  Equation  4.1  then  implies 
that  m1=  m2,  the  classic  Atwood's  machine  case  with  equal  masses.  For  this 

classic  example  there  is  indeed  a  solution  p  =  p0,  a  constant.  If  0  £  0,  then 
Equation  5  implies  that 

2p0  0  =  0 

Equations  3.1  and  6  are  incompatible.  There  is,  therefore,  no  solution  of 
Equations  3  and  4  for  which  p  is  constant. 

We  proceed  now  to  discover  the  relationship  between  m1  and  m2  for  which 
p(f )  and  0(f )  are  periodic  solutions  of  Equations  3  and  4. 

Case  ii)  The  General  Case. 

The  perturbation  technique  requires  that  Equations  3  and  4  be  modified 
and  rewritten  as  follows: 

p  =  e  [m2p02-(m1-m2)  -  !^m202] 

p  0  + 0  =  - 2p0 

where  we  have  introduced  a  parameter  e  and  written  m1=m1/(m1+m2)  and 
m2=m2/(m1+m2).  In  Equations  7  and  8,  we  change  the  time  variable  t  to  a 

pseudo-time  t  by  means  of  the  formula  f=  VuT t  and  get 
p"=  e[m2p0  '2  +  (jj(m2-m1)  -  !^m2a302] 

p0  "  +  idQ  =  -2p  'Q  ' 


where  the  '  indicates  differentiation  with  respect  to  t.  Assuming  that 
p.e.tu  and  are  analytic  functions  of  e: 


P  =  r0*  ♦  €2f2*. .. 

0  =  0Q  +  e01  +  e202  +  ... 

CO  =  (JDq  +  GQJ1  +  G2OJ2  +  .  .  . 

mi  =  M-o  +  6I-L1  +  ^  +  ' * ' 


We  introduce  these  expansions  into  Equations  9  and  10  and  collect  terms 
in  like  powers  of  € .  Since  e  is  an  arbitrary  parameter,  each  coefficient  of 
the  powers  of  e  must  be  zero.  This  results  in  setting  up  an  infinite  system 
of  pairs  of  second  order  differential  equations  which  can  be  solved 
successively  and  recursively.  The  first  set  of  pairs  is  obtained  by  setting 
e=0.  The  two  equations  obtained  are: 

r0"=0  (11) 

roeo  "+  %eo  =  *2roV  <12) 

The  solution  of  Equation  11  is  rQ  =  rQ,  a  constant,  since  we  assume  as  an 
initial  condition  that  r0  '(0)  =  0.  The  solution  of  Equation  12  then  becomes 

0o=0ocos(kt)  (13) 

where  0O  is  the  initial  angular  displacement,  k2  =  u)0/r0  and  we  have  assumed 
that  0O  '(0)  =  0.  It  will  turn  out  that  oj0  plays  no  role  in  the  final  form  of  the 

formulas  for  p  ( t )  and  d(t).  In  terms  of  the  physical  variables,  k2  =  g Ip  (0). 

If  Equation  9  is  differentiated  with  respect  to  the  parameter  e  and  then  e 
is  set  equal  to  zero,  we  obtain: 


ri  '=  m2ro0O  ^0O2) 


(14) 


When  the  values  of  0O  and  0O  ",  as  given  from  Equation  13  and  its  derivative, 
are  substituted  into  Equation  14,  and  remembering  that  r0  =  rQ  ,  a  constant, 
then,  after  some  algebraic  manipulation,  Equation  14  may  be  written  as: 
r1  "=  co0[-|jL0+m2+  (i/4)m20o2-  (3/4)m20o2cos(2kt)] 

Since  we  are  seeking  a  periodic  solution  for  r,  the  secular  term  in  the 
general  solution  of  Equation  15,  the  term  that  would  give  rise  to  a  quadratic 
increase  or  decrease  with  time  in  ,  will  be  eliminated  if  we  set: 

p0  =  m2+  d/4)m20o2 

Otherwise,  r1  will  contain  a  term  of  the  form  (-p0  +  m2+  (i/4)m20o  2  )t2/2 
and,  depending  on  the  relative  magnitude  chosen  for  (jlq  and  (m2+(  i/4)m20o2), 
r1  increases  or  decreses  quadratically  with  time.  The  solution  of  Equation 
15,  for  the  initial  conditions  ^(0)  =  r,{0)  =  0,  is: 


r,=  (3/i6)rom20o2(  cos  (2kt)  - 1) 

When  Equation  10  is  differentiated  with  respect  to  e  and  then  e  is  set 
equal  to  zero,  we  obtain: 


ro0  ,+  0)00,-^0  0-^00-2^00 


Since  rQ,  0Q  and  r,  are  known,  we  substitute  their  respective  values  into  the 
right-hand  side  of  Equation  18  to  obtain: 


rQ0  '+o)o0,=  0o{-[o),  +  (i5/32)m2902  cu0]cos(kt)+ 

[(i5/32)m2eo2co0]cos  (3kt)} 
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(15) 


(16) 


(17) 


(18) 


(19) 
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We  are  interested  in  obtaining  the  particular  solution  of  Equation  19  that 
excludes  the  secular  term.  This  is  done  by  choosing 

uj1  =  -(i5/32)m2e02m0 

With  this  done,  the  particular  solution  of  Equation  20,  with  initial  conditions 
0,(0)  =  0,  '(0)  =  0,  is: 


0,=  (15/256)  0o(m2eo2)[cos  (kt)  -cos  (3kt)] 

Note  that  if  gj,  is  not  chosen  as  in  Equation  20,  then  the  solution  for  0, 

will  go  to  infinity  with  time.  We  interpret  this  by  saying  that  the 
mathematical  procedure  fails.  Thus,  in  order  to  retain  mathematical 
viability,  we  must  pick  u),  as  given  in  Equation  20. 


If  we  were  to  proceed  no  further,  we  would  have  the  following 
approximate  solutions  to  Equations  3  and  4  (after  setting  6  =  1): 

P  =  r0[1  +  (3/i6)m2e02(cos  (2kt)  - 1)] 

0  =  90{cos  (kt)  +  (1 5/256) (m20Q2) [cos  (kt)  -  cos  (3kt)]} 

where  m,=  m2  +  ( i/4)m20o2 

oj  =  (dJI-  (i5/32)m20o2] 
k2  =  iu0/r0 

We  proceed  to  get  the  next  higher  order  terms  in  the  expansions  for 
p,0,m1(  and  oi. 

By  equating  the  coefficients  of  e2from  both  sides  of  Equation  9,  we  get 

r2  -[oj0h,  +  uj.qj.,  -  m2)]  +  m2[r,e0'?  +  pr0e0'e,'] 


After  substituiting  the  several  quantities  already  found  for  their  respective 
values  as  given  in  the  right-hand  side  of  Equation  27,  it  may  be  rewritten  as: 

r2  "  =  (9/512)m220o4Q3o  [4  cos  (2kt)  +  9  cos  (4kt)  ] 

where,  in  order  to  preclude  the  introduction  of  a  secular  term  in  r2  ,  we  had 
set: 

p,  =  (63/51 2)m22e01 

The  solution  of  Equation  28,  subject  to  the  initial  conditions  r2(0)  =  r2  (0)  = 
0,  is: 

r2  =  (9/213)(m20o?')2ro[25  - 16  cos  (2kt)  -  9  cos  (4kt)] 

By  equating  the  coefficients  of  €2from  both  sides  of  Equation  10,  we  get: 

roe  (jJoe?  +  [r,©  ''  +  (Dj ©j  +  2r1  ']  + 

[r2e'0'+  iu2e0  +  2r2  'ea']  -o 

After  substituting  the  values  already  found  for  ro,0o,r1 ,01  ,oo0 .oj,  and  their 
derivatives  where  appropriate,  we  get: 

ro0,. '  +  a)o02  =  (9/214)m22©05cu0[i63  cos  (3kt)  -  291  cos  (  5kt)  ] 

and  where,  in  order  to  preclude  the  introduction  of  a  secular  term  in  02,  we 
have  set: 

uj2  =  (9/i28)m220o4coo 

The  value  of  02  ,  subject  to  the  initial  conditions  02(O)  =  0^(0)  =  0,  found  by 
integrating  Equation  (32)  is: 


02  =  (9/217)m22  625  [  66  COS  (kt)  -163  COS  (3kt)  +  97  COS  (5kt)] 

Summarizing  the  results  obtained  thus  far,  we  see  that  after  setting  e  = 
1 ,  we  get  the  following  approximate  solutions  for  p(t)  and  0(t): 

p(t)  =  r0Ji  +  (3/i6)m2e22(cos  (2kt)  -i)  + 

(9/213)m220o4(25-16COS  (2kt)-9  COS  (4kt)] 

0(t)  =  0o{cos  (kt)  +  (l5/256)m20Q2[  cos  (kt)  -  cos  (3kt)]  + 
(9/217)m22024  [  66  cos  (kt)  -163  cos  (3kt)  +  97  cos  (5kt)]}, 

where 

m1  =  m2+  (i/4)m20Q2  +  (63/51 2)m220o4 
m  =  o)0[i-  (i5/32)m20o2  +  (9/i28)m220o4] 

Since  t  =  (id) m*t,  p  (t)  =*  p  (0) )  ,  0(t)  =>  0  (oT^  t )  ,  i.e.  in  formulas 
(35)  and  (36)  replace  t  by  uT^f.  Finally, 

p(f)  =  gp(ai^f) 

and 

Q(t)  =  0(w  ^f) 


Section  III 


Discussion 


An  understanding  of  the  results  presented  in  the  previous  section  may  be 
enhanced  by  a  numerical  example.  Before  doing  this,  however,  we  want  to 
recall  the  relationships: 

p(f)  =  gp(f)  and  f  =  vrurt 

These  imply  that  p  =  g  p  =  g  p"/w  .  Using  these  we  can  calculate  the 
acceleration  acting  on  the  mass  m^ .  We  do  this  first  using  the  approximate 

solutions  given  by  Equations  22  and  23.  In  particular,  we  examine  Equation 
22  and  take  the  second  derivative  with  respect  to  t  to  get: 

p"(t)  =  -(3/4)m2e02(D0cos  (2kt)  (39) 

and  where  we  have  used  k2  =  cjo0/r0  .  Then,  using  the  value  of  co  as  given 
in  Equation  25  : 


p  =  -{(3/4)[  m2e02]/  [1-  (i5/32)m2602]  }cos(2kt)  (40) 

Remembering  that  m1=  1  -  m2,  we  solve  Equation  24  for  m2  and  substitute 
that  value  into  Equation  (40)  to  get: 

p  =  -[240o?/(64  -  70*)]  cos(2kt)  (41) 

For  values  of  0O?  of  interest,  say  not  to  exceed  0.1  radians,  the  expression 
within  the  [...]  in  Equation  41  is  a  monotonically  increasing  function  of  @0  and 
will  therefore  attain  its  maximum  value  at  the  end-point  of  the  interval  of 


(42) 


interest,  if  we  take  that  to  be  0  <  0O  <  0.1  radians,  then  for  60  =  0.1 , 

'p(t)  =  -.0037541  g  cos  (2kt) 

Thus,  the  acceleration  of  the  mass  m1  is  only,  at  most,  0.00375  times  the 
acceleration  of  gravity,  a  very  small  quantity.  While  the  mass  m1  does  not 

move  as  if  it  were  in  free  space,  it  is  subject  to  a  gravitational  influence 
that  is  only  .00375  times  the  effect  it  would  normally  sense  on  earth. 

The  description  given  above  is  based  on  using  only  two  terms  in  the 
expansion  for  p  (t ).  Some  change  in  the  calculated  value  of  the  acceleration 
of  the  mass  m1  should  be  expected  if  three  terms  are  used  in  the  expansion 

for  p  ( t ).  While  it  is  possible  to  carry  out  these  computations 
algebraically,  no  new  insight  into  the  motion  is  gained.  Accordingly,  so  as 
not  to  get  lost  in  a  maze  of  algebraic  manipulation  and  symbolism,  we  limit 
our  discussion  to  a  numerical  example,  indicating  along  the  way  how  other 
numerical  examples  may  be  handled. 

We  start  with  Equation  35  and  calculate: 

p"(t)  =  qj0  (m20o2)  [A1  cos  (2kt)  +  A2  cos  (4kt)]  (43) 

where: 

A1  =  [-(3/4)  +  (9/128)  n^Gg2] 

A2  =  (81/512)  m20o2 

p 

and  again  we  have  used  k  =ajQ/r0  .  Then: 

p  (t)  =  A0[A1  cos  (2kt)  +  A2  cos  (4kt)]  (44) 

where: 

A0  =  (m2e02  uj0)/  cjj 

or: 
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A0  =  m2e02/[i-(i5/32)m2e02  +(9/i28)m22e04] 


$ 


k 

a 


R 


where  the  value  of  (D  is  taken  from  Equation  38  . 

The  acceleration  is  stationary  when  2kt1=0+nTr  and  when  cos  (2kt2)  = 
-A1/(4A2).  For  these  values  of  t1 ,  the  maximum  value  of  the  acceleration 
cannot  exceed  |A0|g[  |AJ  +  |A2|] ;  for  values  of  0O  of  interest,  the  magnitude 
of  A/(4A2)  is  greater  than  1  and  there  is  no  real  value  for  t2. 

In  order  to  calculate  values  for  A0,  Av  and  A2  ,  we  need  to  know  m20o2. 

For  a  given  value  of  0Q  ,  say  0Q  =  0.1  radians,  and  remembering  that  m1  +  m2 
=1 ,  we  find,  using  Equation  37,  that  m2  =  .4993742479515557.  We  then  find 
that  A0  =  5.0054505  x  10'3,  A1  =  -.74964888,  A2  =  7.9002567  x  10'4  and  the 
maximum  value  of  the  acceleration  is  at  most  3.7483759  x  10'3g. 

A  comparison  of  the  results  when  2  or  3  terms  are  used  in  the 
determination  of  the  acceleration  of  m1  indicates  that  the  difference  is  so 

small,  of  the  order  of  5.7  x  1 0‘ 6  g,  as  to  be  negligible.  For  all  practical 
purposes,  the  solution  to  the  problem  is  given  by  the  Equations  22  -  26  and 
Equation  41 . 

When  the  fundamental  frequency  of  the  system  is  compared  to  the 
frequency  of  a  pendulum  of  length  p  (0),  there  is  a  change  in  frequency  given 
by  the  multiplicative  factor  [(64  +  80o2)/(64  -  70Q2)]1/2.  For  ©0  =  0.1  radians, 
this  factor  is  1 .001 1 725. 


The  solution  to  Equations  1  and  2  as  given  by  Equations  22-26  or,  if  you 
prefer,  by  Equations  35-38,  are  periodic  solutions.  As  already  noted,  to 
obtain  a  periodic  solution  for  p  it  was  necessary  to  carefully  adjust  the 


ratio  of  the  two  masses  m1  and  m2. 


If  values  of  m1  and  m2  different  from 


those  determined  by  means  of  Equation  22  (or  Equation  37)  are  used,  then  the 
solution  for  the  pendulun  arm  p  (f )  will  contain  an  additional  term  involving 


t 2 .  The  case  for  m1  =  m2  was  studied  in  reference  [8]. 

If  a  value  of  w  different  from  the  ones  given  in  Equations  25  or  38  is  used, 
then  the  mathematical  process  fails  because  secular  terms,  which  physically 
and  experimentally  we  know  do  not  have  meaning,  are  introduced  into  the 
equation  for  d  ,  Equations  23  or  36. 
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Section  IV 


Conclusions 


A  theoretical  description  of  an  Atwood's  pendulum  machine  is  given. 

Implicit  in  these  calculations  were  the  assumptions  that  the  pulley  is 
frictionless  and  that  its  radius  could  be  neglected.  Also,  the  wire  connecting 
the  two  masses  is  assumed  to  be  perfectly  flexible,  inextensible  and  of 
negligible  mass.  The  following  results  were  deduced: 

(1)  If  the  ratio  of  the  two  masses  are  chosen  as  prescribed,  then  the 
vertical  motion  of  the  mass  m1  ,  which  is  the  only  motion  possible  for  this 

mass,  is  periodic; 

(2)  Equivalently,  the  change  in  length  of  the  pendulum  arm  associated 
with  the  swinging  of  the  mass  m2  is  exactly  the  same  as  the  change  in  length 

of  the  wire  supporting  m1  except  that  as  the  one  increases  the  other 
decreases; 

(3)  If  the  two  masses  are  not  picked  as  prescribed  so  as  to  insure 
periodic  motion,  then  the  motion  of  m1  consists  of  the  superimposition  of 

two  motions: 

(a)  an  accelerated  motion  as  if  the  mass  m1  were  in  a 

gravitational  field  of  magnitude  <  g,  the  magnitude  of  this  deviation  from  g 
depending  on  just  how  much  the  two  masses  differ  from  the  ratio  required 
for  the  periodic  motion,  and 

(b)  the  periodic  motion  superimposed  (just  added  onto)  the  motion 
described  in  (a); 

(4)  While  the  conclusions  (1),(2),  and  (3)  given  above  were  obtained  by 
considering  only  two  terms  in  the  expansion  of  the  relevant  functions,  the 
effect  of  using  three  terms  in  these  expansions  was  shown  to  modify  the 
conclusions  only  quantitatively,  not  qualitatively,  and  the  change  was  of  an 
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In  mathematical  terms,  this  periodic  solution  is  referred  to  as  a  limit 
cycle;  it  is,  however,  unstable.  Any  set  of  conditions  that  differ  ever  so 
slightly  from  the  ones  needed  to  establish  the  limit  cycle  will  produce  a 
trajectory  that  increasingly  diverges  from  the  limit  cycle  and  this 
translates  into  the  statement  that  the  pendulum  arm  goes  either  to  zero  or 
to  infinity^10'. 


Figure  1 

A  Schematic  for  an  Atwood's  Pendulum 
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List  of  Symbols 


=  mass  1 


=  mass  2 


=  reduced  massl  =  m1/(m1+m2) 

,  =  reduced  mass  2  =  m^rr^+m^ 

=  real  time,  seconds 
=  pseudo-time 

=  initial  angular  displacement 
=  initial  angular  velocity 
=  length  of  pendulum  arm,  meters 

=  6  =  angular  displacement  of  pendulum  arm;  two  symbols  are  used  to 
make  the  notation  uniform 

=  p/g  =  length  of  pendulum  arm  in  a  system  of  units  in  which  g  =  1 
=  a  parameter  relating  t  and  t;  eventually  shown  to  be  a  factor 
specifying  frequency 

=  a  parameter  used  to  establish  asymptotic  expansion 
=  a  term  in  the  expansion  of  p 

=  a  term  in  the  expansion  of  0 
=  a  term  in  the  expansion  of  w 
=  a  term  in  the  expansion  of  m1 

=  initial  length  of  pendulum  arm  in  system  of  units  for  which  g  =  1 


a  parameter  relating  to  frequency  in  the  system  in  which  g  =  1 ; 

,,,  / , 
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